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Based on the stochastic superspace mechanism for softly breaking supersymmetry, we present 
improved sparticle spectra computations for the minimal model and examine extensions through 
R-parity violation and the type-I seesaw mechanism that incorporate non-zero neutrino masses 
for more realistic models. Performing the calculations to two-loop accuracy, we observe a global 
decrease in predicted sparticle masses. However this does not affect the generic features of the 
minimal model outlined in our earlier work, including the characteristic light stop mass. We find 
stop decay channels accessible at the LHC which can be used in combination with our predicted 
range for the stop mixing angle to falsify the minimal model with stochastic supersymmetry. We 
then introduce neutrino masses and mixings consistent with experiment by including purely trilinear 
R-parity violating superpotential terms, resulting in a viable stochastic superspace model absent a 
dark matter candidate. An alternative method for generating neutrino masses, namely the type-I 
seesaw mechanism, is found only to be viable when the neutrino Yukawa coupling is small relative 
to the top Yukawa and the cut-off scale is large. 
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I. INTRODUCTION 



A number of methods exist for softly breaking supersymmetry. These can be broadly catagorized into two groups: 
Oh' those that originate from a fundamental theory, where supersymmetry is spontaneously broken in a hidden sector then 
communicated to the visible sector by a messenger, and the purely phenomenological approach, where no explanation 
• for the cause of supersymmetry breaking is offered. 

The fundamental theoretical approaches, though providing some dynamical explanation for the mechanism of 
l— ~~ '■ supersymmetry breaking manifesting in the observable sector, suffer from a variety problems. Gravity mediation 
is plagued by large sparticle induced FCNC amplitudes unless no source of flavor physics exists between the weak 
scale and Planck scale, while gauge mediation has CP violation issues among other significant theoretical hurdles 
£SJ ' Although some work-arounds for these problems exist, there is not yet a compelling and simple fundamental 

model for the transmission of supersymmetry breaking to the observable sector. On the other hand, writing down 
the most general Lagrangian with explicit soft supersymmetry breaking terms completely fails to explain the source 
of supersymmetry breaking. Instead, one must painstakingly explore the vast parameter space, much of which is 
£f~) ■ excluded due to flavor changing neutral current (FCNC) and CP violation problems, to find phenomcnologically 
viable regions useful for experimental predictions. 

We proposed in [2j a mechanism for supersymmetry breaking in which the Grassmannian coordinates, 9 and 9, 
are considered to be stochastic variables. This leads to a very constrained set of soft-breaking parameters, thus 
significantly improving predictability compared to the phenomenological approach. However, the underlying cause of 
stochasticity in the Grassmannian coordinates is not postulated, placing stochastic superspace as a sort of compromise 
between a fully dynamical model and the purely phenomenological approach to supersymmetry breaking. 

One of the goals of this paper is to improve upon the sparticle spectra calculations of the minimal stochastic 
superspace model presented in [2] . Rather than using the analytic solutions to approximations of the one-loop renor- 
malization group equations (RGEs), we employ the use of the publicly available software, SDFTSUSY3.0 [3|, to account 
for two- loop and threshold effects. We then incorporate non-zero neutrino masses into the stochastic superspace mech- 
anism for supersymmetry breaking in order to create more realistic models. We examine the phenomenology of models 
that achieve this end through R-parity violation and the type-I seesaw mechanism. 

The rest of the paper is structured as follows: we briefly review the stochastic superspace formalism in the next 
section and present updated sparticle spectra. We then apply the stochastic superspace formalism to R-parity violation 
(RPV) and type-I seesaw in order to generate non-zero neutrino masses in sections IIIII and IIV1 respectively. This is 
followed by concluding remarks in section IVl 
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II. THE STOCHASTIC SUPERSPACE FORMALISM 

In stochastic superspace models the Grassmannian coordinates, 9 and 9, are considered to be stochastic variables 
with a probability distribution, 

V(9, 0)|£| 2 = 1 + e{66) + £(00) + |£| 2 (00)(00), (2.1) 

which is dependent on a single parameter, £, a complex number of mass dimension. This probability distribution 
is the most general form that satisfies Lorentz invariance. The Lagrangian in ordinary spacetime is then defined 
as the average of the supersymmetric Lagrangian over the probability distribution in Eq. (|2.1j) . Averaging over the 
superpotential for the minimal supersymmetric standard model (MSSM), 

Wsm = ^H u H d + y up QU c H u + y down QD c H d + y lcpt LE c H d , (2.2) 

gives rise to the soft-breaking terms 



isoft scalar — — C^H u H d — 2£* 



r P QU c H u +y dom QD c H d + y lopt LE c H d +H.c, (2.3) 



with the tildes denoting the scalar component of the chiral superfield in question. Similarly, averaging over the the 
gauge-kinetic F densities lead to soft-breaking masses for the gauginos, X^(x), 



igaugc 



i^TrVF (l)Q W^ : 



- ^^TrA (l) A (l) +H.c, (2.4) 



where W^ a represent the standard model field-strength superfields. Refer to @ for further details of the derivation. 

At some cut-off scale, A - also a parameter of the theory - the soft-breaking terms for the minimal stochastic 
superspace model are: 

B» = e, (2.5a) 
A = 2C, (2.5b) 

mi/2 = ifcl, (2.5c) 

ml = 0. (2.5d) 

Here, is the bilinear Higgs soft term, Aq is the trilinear universal coupling, toi/2 is the universal gaugino mass 
and mo denotes the universal soft scalar mass. Thus, the minimal model in stochastic superspace represents a further 
constrained version of the so-called constrained MSSM (CMSSM). The analysis of this model can be summarized by 
stating the key phenomenological features at the extrema of acceptable cut-off scales: 

• For A = Mqut, the lightest supersymmetric particle (LSP) is the stau. Since this model conserves R-parity, 
the LSP is a stable particle. Thus, this model is excluded by experiment. 

• For A = Mp\, the LSP is the neutralino. A phenomenologically viable region of ^-parameter space exists at this 
cut-off scale. 

However, as this model of soft supersymmetry breaking from stochastic superspace is based on the MSSM, neutrino 
masses are inherently absent. In the following sections we explore the possibility of applying the stochastic superspace 
formalism to two of the standard mechanisms for generating neutrino masses in supersymmetric models: through R- 
parity violation and type-I seesaw. While doing so, we aim to find phenomenologically viable regions of parameter 
space that do not necessitate pushing the cut-off as high as the Planck scale, where quantum gravitational effects 
become significant. 



A. Improvement of the Sparticle Spectrum Calculation 

In the original analysis of the stochastic superspace mechanism for soft supersymmetry breaking, the calculation 
of the superparticle spectrum was performed using an analytical solution to an approximation of the one-loop renor- 
malization group equations 0, 0|. We have since improved upon the accuracy of these results by performing the 
calculations using the sparticle spectrum calculator software SDFTSUSY3.0 [1, @, which uses an iterative numerical 
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technique to solve the two- loop RGEs, whilst taking into account threshold effects. Figure [T] displays sparticle spectra 



calculated by S0FTSUSY3.0 for various points in the (£, A) parameter space, with Figure 1(b) corresponding to the 
same parameter choice as the example spectrum presented in Q . Of note is a global decrease in the predicted sparticle 
masses, but this does not affect the generic features of the minimal stochastic superspace outlined earlier. For example, 
where the neutralino is LSP, there is still a large region of parameter space for which the lightest stau is no heavier 
than 10 GeV more than the neutralino mass. This is the stau coannihilation regime, which ensures that processes in 
the early universe result in the correct cosmological dark matter abundance. Additionally, the character of the LSP 
swaps between stau and neutralino when A ~ 6 x 10 17 GeV, depending slightly on £. x This is also consistent with 
our earlier results. 

As has been noted earlier, the minimal model with stochastic supersymmetry is a particular case of the CMSSM. 
Therefore, besides sharing generic phenomenological features with the CMSSM, our model gives more specific predic- 
tions. Namely, tan/3 is no longer a free parameter, but predicted to be tan/3 « 4-^5 (depending on £ and A). Despite 
such a low tan/3, the mass of the lightest CP-even Higgs boson is predicted to be ps 112 115 GeV due to the 
enhanced loop contribution stemming from the large mixing parameter in the stop sector. 2 

Another characteristic feature of the stochastic superspace sparticle spectra is a lighter-than-usual stop mass. 3 
When the cut-off scale is close to Mpi, the stop decays predominantly through the process, 

k -> tx°v (2.6) 

If, say, the LHC were to detect a larger cross-section for top production than expected, it could suggest stochastic 
superspace with a high cut-off scale. While this is true for £ < —500 GeV, as £ becomes more positive the stop decays 
primarily through the also experimentally interesting channel, 

h -> bxf. (2.7) 

The stochastic superspace hypothesis (at least the minimal model) can be falsified if none of these signatures are 
observed in future experiments. However, if (|2.6|) is observed at LHC, one might be able to extract the important 
stop mixing angle 6* from the measurements of forward-backward asymmetries in leptonic and hadronic decays of 
polarized top quark 0. Also, simultaneous analyses of decays (|2.6p and (|2.7p will allow 9t to be determined with high 
accuracy [lfj. Our model predicts: 9t ~ 0.3 -j- 0.4. 4 

Finally, we have analyzed the amount of fine-tuning in our model by computing 



A = max 



<91ogm| , ^ 



d log p t 



where pi is a relevant parameter [11| . The fine-tuning measure is found in the range A w 200 -j- 900, with the lower 
limit corresponding to £ = —500 GeV and the upper limit corresponding to £ = —1000 GeV. This means that the 
fine-tuning is 2% -f- 10% over the allowed (£, A) parameter space. 



III. R-PARITY VIOLATION 



In supersymmetric extensions of the standard model it is generally possible to find renormalizable, gauge-invariant 
interactions that violate either lepton number or baryon number conservation. This is particularly problematic when 
superpotential terms that violate lepton and baryon number conservation are simultaneously present, as this will lead 
to disastrously rapid decay of the proton. The MSSM imposes a further symmetry, called R-parity [l2|, [l3[ , to forbid 
such dangerous phenomenological effects from arising. R-parity is a discrete symmetry defined by 

R p = (_1)2S+3B+L ; (3.!) 



1 In the literature the region of parameters with mo = is usually considered to be excluded experimentally, due to the fact that relatively 
light stau is the LSP. This is simply because the soft- breaking parameters are customarily defined at A = A/gut- Ellis et al. have also 
recently noticed in [(J that this region of parameters is acceptable with higher cut-off scales in the CMSSM, consistent with our earlier 
results in Q|. 

2 Taking into account typical 3-^5 GeV theoretical error in determining m/, Q , our prediction is consistent with the LEP II lower bound 
on SM-like Higgs boson, m^ EP > 114.4 GeV 

3 Light stops are welcome within the context of electroweak baryogenesis, since they trigger a strong first-order electroweak phase transition 

4 \0t\ < tt/4 is defined as in [3. 
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(e)A = Mqut.C = -1000 GeV. 
FIG. 1: Sparticle spectra computed using S0FTSUSY3.0, for various values of the parameters A and £. 



so that standard model particles are R-parity even (R p = +1) and their superpartners are R-parity odd (R p = — 1). 
When R-parity is conserved, R-parity odd particles can only be produced in pairs, and sparticle decays must produce 
an odd number of sparticles. The phenomenological benefits of R-parity conservation are thus two-fold: lepton number 
and baryon number are conserved, and stability of the lightest supersymmetric particle (LSP) is guaranteed, which 
often provides a natural dark matter candidate. 

From a theoretical standpoint, however, the choice of R-parity conservation is arbitrary. If we permit R-parity to 
be violated, the superpotential should be augmented by the terms, 



- e af3 ~< A " jk Ul a Djp D c kl , 



(3.2) 
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where {to, n} are SU (2) l indices and {a, [3, 7} are SU (3) gauge indices. The first three terms violate lepton number 
conservation, and originate from replacing Hi in the R-parity conserving superpotenial with Li, which has identical 
gauge quantum numbers. The last term violates baryon number, and is the only remaining renormalizable and 
gauge-invariant superpotential term allowed when considering the same field content and gauge symmetries as the 
MSSM. 

It is still possible to avoid inconvenient effects such as the aforementioned proton decay if a less constraining 
symmetry - baryon triality - is considered [14| . Here, the last term in Eq. (|3.2[) is forbidden, effectively allowing 
only lepton number violation. One of the key phenomenological differences resulting from R-parity violation is the 
generation of neutrino masses and mixings. Additionally, since the LSP can decay into standard model particles, its 
charge or colour neutrality is irrelevant because it no longer presents as a dark matter candidate. 



A. RPV and Neutrino Masses 



The presence of bilinear R-parity violating terms in the superpotential leads to an interaction between leptons and 
higgsinos that was otherwise absent in the MSSM. These new interactions require the neutralino mass matrix to be 
extended to a 7 x 7 matrix, which includes the three neutrinos in addition to the four neutralinos [l5j . The structure 
is reminiscent of the seesaw mass matrix, and the same techniques can be used to extract the 3x3 neutrino mass 
matrix. The tree-level neutrino mass matrix, 

-.2 



M cos 2 P 

[""Itrcc-levelJi,- ~ 77 MiMj: IA<*) 

M SVSY 



results in a single non-zero eigenvalue. Mgysy denotes the scale at which supersymmetry is broken. Solar and 
atmospheric data for neutrino mass-squared differences [16| indicate that at least two neutrinos should be massive. This 
rules out the possibility of purely bilinear RPV interactions generating neutrino masses consistent with observations. 
Trilincar RPV superpotential terms couple a scalar field with two fermionic fields, resulting in radiatively generated 
neutrino masses from slepton- lepton and squark-quark loops. For purely trilinear effects 5 , the loop contributions can 
be expressed as (l7| . 

r ,(AA) 1 \ \ m lk m h to A\ 

ott jWsusy 



[mioopnj 



k.i 



lir 2 lKl ^ AfsUSY 



where mi k and rrid k are the fc th generation charged-lepton and down-quark masses, respectively. The neutrino mass 
matrix is then the combination of the bilinear tree-level contributions with the trilinear loop corrections, 

[m v ]ij = [TO troe _ievei]-^ M) + [mioop]y + [mioop]^ X \ (3.6) 

which can, in general, yield three non-zero eigenvalues. In order to find an R-parity violating model that produces 
neutrino masses and mixings consistent with experimental data, we choose to diagonalize the mass matrix in Eq. p.6[) 
with the tribimaximal mixing matrix, defined by 




1= l-Jj, ^ (8.7) 
75 

using the phase convention of (T^ |. The neutrino mass matrix can then be written in terms of the mass eigenvalues, 
to„ = J7diag(mi, nx%, m^)U T 

-TO1+TO2) i(-TOl+TO 2 ) \ 

2to 2 + 3m 3 ) I (toi + 2to 2 - 3to 3 ) . (3.8) 

2to 2 - 3m 3 ) g (to! + 2to 2 + 3TO 3 )y 




5 The bilinear RPV interactions can also generate loop contributions which may dominate under certain conditions. See [17| for a summary 
of other possibilities. 
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Parameter 


Normal Hierarchy 


Inverted Hierarchy 


A133 


8.058 x 10" 5 


- + - 


+ 


3.503 x 10" 4 


1 h 


A233 


1.612 x 10~ 4 


- + - 


+ 







A'l33 


6.802 x 10~ 6 


h 


+ 


6.339 x 10" 7 


h + 


^233 


5.035 x 10~ 5 


+ + - 




6.077 x 10 -5 


h + 


-^333 


6.395 x 10 -5 


h 


+ 


6.077 x 10" 5 


h + 



TABLE I: R-parity violating parameter sets (at weak scale) that achieve tribimaximal mixing in the normal hierarchy (mi = 0, 
m 2 = 8.75 x 10" 12 GeV, m 3 = 4.90 x 10" 11 GeV) and the inverted hierarchy (mi = 4.90 x 10" 11 GeV, m 2 = 4.98 x 10" 11 GeV, 
?7i3 = 0). The columns of '+' and ' — ' represent the eight possible cases, indicating the sign of the corresponding parameter 
whose magnitude does not change for a given choice of mass hierarchy. 



Setting this equal to the expression in Eq. (|3.6l) defines sets of R-parity violating couplings at the weak scale that give 
rise to experimentally consistent neutrino masses and mixings, provided an appropriate scale for the neutrino masses 
is chosen. 

As our goal is to establish whether or not a phenomenologically permissible realization of R-parity violation is 
possible in the stochastic superspace framework, we proceed to analyze a constrained sub-section of the vast parameter 
space offered by the RPV regime. Contingent on a positive result, we argue that a thorough scan of a more relaxed 
parameter space could map further regions compatible with stochastic superspace. However, a study of this magnitude 
is beyond the scope of this work. 

For simplicity, we set one neutrino mass (which fixes the value of the remaining two from solar and atmospheric 
data) and all bilinear couplings at weak scale to zero 6 . Inspection of the loop contributions in Eq. (|3.4l) and Eq. p.5[) 
reveals a dependence on products of lepton or down-quark masses, respectively. Given the mass hierarchies in these 
sectors, we make the approximations m T ^> m^, m e and mi, ^ m s , mj, which simplifies the loop contributions to 



[m loop ]| AA) » 1 " v Ai3 3 Aj33 (3.9) 

87T Z MsuSY 



and 



[-loopllf } « j^T^-WW ( 3 - 10 ) 

Effectively, the relevant trilinear R-parity violating parameter set has been reduced to 7 

{A133, A233, A 133 , A 23 3, A333}. (3-11) 

The five resulting independent linear equations can then be solved for the set of weak scale trilinear RPV couplings 
that achieve neutrino masses in the normal or inverted hierarchy with tribimaximal mixing. See Table U for the eight 
possible parameter sets for the case where one neutrino mass is zero. 



B. Results 



Determining the phenomenology of soft supersymmetry breaking from stochastic superspace in the context of 
R-parity violation follows the same procedure as for the R-parity conserving case: set the renormalization group 
equations' boundary conditions at some scale A using the conditions in Equations (|2.5[) and run down to the weak 
scale to obtain the sparticle spectrum. We encounter a complication in the R-parity violating case in that we require 
a very specific set of RPV trilinear couplings at the weak scale to obtain realistic neutrino masses and mixings. These 
weak scale values must be translated to a set of cut-off scale inputs. 

We find empirically, using the RPV sparticle spectrum calculator software S0FTSUSY3.0 Q, that there is an ap- 
proximately linear relationship between the weak and cut-off scale values of the trilinear couplings at the order of 



Due to the nature of the RGEs for fii, fii(M z ) = does not imply the same condition holds at some higher cut-off scale. Refer to the 

review [19] for a detailed explanation. 

A333 = from antisymmetry in the first two indices. 
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magnitude of interest for neutrino mass generation. The proportionality constant differs for charged-lepton and down- 
quark couplings, and is also affected by the mSUGRA cut-off scale inputs. It must therefore be recalculated at each 
point in the (£, A) parameter space. We define, 

\{M Z ) ee 6A(A), (3.12) 
\'{M Z ) ee b'A'(A), (3.13) 

where the constants b and b' are found by setting all but one of the RPV couplings to zero. We do not discriminate 
between individual couplings within the charged-lepton and down-quark sectors since the difference is negligibly small. 
It is then straightforward to find the cut-off scale input values for the trilinear R-parity violating couplings. Non-zero 
trilinear couplings, however, generate contributions to the bilinear RPV couplings through the running of the RGEs. 
The contributions from the lepton and down-quark couplings are approximately independent, so we can express the 
bilinear coupling at weak scale as a linear combination of the trilinear contributions, 

m(M,) = aA l33 (A) + a% 3 (A), (3.14) 

where the constants a and a' are defined by the ratios, 



A(A) 



(3.15) 



a = (3.16) 

This would result in non-zero bilinear RPV couplings at the weak scale, which would then affect our previous analysis 
of the neutrino masses, where we have neglected the contributions from these couplings by assuming they are non- 
existent at the weak scale. In the absence of trilinear couplings, we notice that, like the trilinear couplings, there is 
an approximately linear relationship between the weak and cut-off scale values such that, 

/*(Af,) ee d/j,(A). (3.17) 

It is conceivable, then, to engineer a situation where the bilinear couplings at the cut-off scale have been chosen so 
that the contribution from the trilinear couplings in Eq. (|3.14[) is exactly cancelled at the weak scale, i.e., we set 

Mi(A) = _aA i3 3(A) + a^ 33 (A) (3ig) 

Using a modified version of S0FTSUSY3 .0.11 and the above algorithm for determining the appropriate RPV cou- 
plings at the cut-off scale, we scanned the stochastic superspace parameter space in the region —1000 < £ < —100 
and Mqut < A < Mpi. Because of the small magnitude of the RPV couplings, their contributions to sparticle 
masses through new terms in the RGEs 8 are negligible, so the sparticle spectra for R p conserving and R p violating 
stochastic superspace models are not distinguishable (see Figure [l}. Like in the minimal stochastic superspace model, 
we find that the stau is LSP for A = Mqut, and as the cut-off scale is pushed higher towards the Planck scale, it 
changes to the neutralino; see Figure [5] The crucial difference is that in the R-parity violating scenario the LSP is no 
longer stable, so a stau LSP is not automatically excluded by experiment. It decays into Standard Model particles 
predominantly through the two-body processes [20(, 

Ti^T-Vi, (3.19) 
Tf-tT-Vi, (3.20) 
ff -> lrv T . (3.21) 

The branching ratio for each of these processes is equal, with the 2-body decay width given by 

r \ - 3|A i33 | 2 mf . ooa 

1 two-body ~ ' (o.ZZ) 



8 See the appendix of [T^] for a full listing of the RPV RGEs used in the SDFTSUSY3 . software. 
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(a)A = M G ut- For £ < -300 GeV, the stau is the LSP. (b)A = M P1 . The neutralino is the LSP. 



FIG. 2: These plots display the £— dependence of the masses of the lightest stau (m T i), lightest neutralino (m x i), lightest stop 
(mti) and lightest CP-even Higgs for two representative choices of the cut-off scale, A. 



which puts an upper bound on the lifetime of the stau LSP of the order 10~ 18 s, guaranteeing its decay occurs inside 
the detector. Contrasting to the decay channels in Eqs. (12.6^ - (|2.7|) . when R-parity is violated the stop can decay 
directly into Standard Model particles through the couplings A- 33 . The dominating process, 

h -)• br, (3.23) 

should be within experimental capability for detection at the LHC. 

Flavor-changing processes are often an issue with R-parity violating models. However, due to the specific flavor 
structure of the R-parity violating couplings detailed in Eq. p. lip , interesting flavor-changing processes are not induced 
at tree-level. Many are induced indirectly through the generation of relevant couplings due to the RG running. 
However, the couplings generated in this way are smaller than the ones in Tableland thus they easily satisfy existing 
experimental bounds. There are also some radiative flavor-changing processes which are directly generated in our 
model. For example, in the model with normal neutrino mass hierarchy, the dominant contribution to the radiative 
decay /i — > comes from the A couplings. The branching ratio is estimated as: 



Here a w 1/137 is the fine-structure constant and Gf ~ 1.166 • 10 -5 GeV -2 is the Fermi constant. In the case of 
inverted hierarchy, A couplings do not contribute directly and the dominant contribution is given by A' couplings. We 
find the branching ratio to be: 

B^e 7 )*|^^~10-x(^y. ,3.25, 
16^(G f t%) 2 V m i n J 

In both cases the predicted branching ratios are many orders of magnitude smaller than, not only the current ex- 
perimental upper limits, but also upper limits on the branching ratio (~ 10~ 14 ) which can be reached in future 
experiments. A similar situation pertains to other flavour-changing processes such as /i — » eee decay and fi — > e 
conversion in nuclei. 

We conclude that the R-parity violating extension to the minimal stochastic superspace model is phenomenologically 
viable, producing neutrino masses and mixings consistent with experiment whilst keeping the magnitude of the RPV 
couplings small enough that they do not conflict with published constraints [Til [2l| . The only outstanding issue is 
the removal of a dark matter candidate from the theory. This can be remedied by, for example, the introduction of 
axions and axinos. We do not address this here, but aim to explore the details of RPV dark matter with stochastic 
superspace in future work. 



IV. TYPE-I SEESAW 



Another well-motivated way to introduce neutrino masses to a supersymmetric model is through the seesaw mech- 
anism, of which there are three types. We base the following analysis on the simplest of the three: the type-I seesaw 
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mechanism. For type-I seesaw, the held content of the R-parity conserving MSSM is extended by three additional 
gauge singlet chiral superfields, N, for the right-handed neutrinos which we assume to be degenerate. The seesaw 
superpotential, 

= f iH u H d + y np QU c H u + y dom QD c H d + y lcpt LE c H d + y ncut LN c H u + hl R N c N c , (4.1) 

includes a neutrino Yukawa coupling and Majorana mass term that were not present in the MSSM. A neutrino Dirac 
mass matrix arises from the new Yukawa term, giving 

m D = y mut (H u ), (4.2) 

where (H u ) — v sin/3 is the vacuum expectation value for the Higgs held with positive hypercharge, v = 174 GeV and 
we make the dominant third family approximation with respect to the Yukawa couplings. The light neutrino mass 
matrix, 

m v — —mj^M^niD, (4.3) 

follows from evaluating the eigenvalues of the neutrino mass matrix under the assumption that the Majorana mass 
scale is much greater than the Dirac mass. Adopting the notation of 3, we define the largest coupling 



fw ncut r 

l07T^ 



which becomes a free parameter of the theory. The magnitude of Y/v sets the seesaw scale, 

M R ~Y n x 10 17 GeV, (4.5) 

where we have assumed the light neutrinos' masses to be O(10 _11 GcV). For simplicity, in the following section we 
strive to generate light neutrinos with mass of this scale, without taking into account the experimentally observed 
mixing and mass difference data. 



A. Results 



Introducing the right-handed neutrino superfields affects the RGEs above the seesaw scale; see Appendix [5] for a 
listing of the equations that differ from the MSSM. To find the sparticle spectrum for our model with soft breaking 
from stochastic superspace, we set the cut-off scale boundary conditions as per Equations (|2.5I) and run the equations 
down to the seesaw scale, Mr using the modified RGEs. Below the seesaw scale, the right-handed neutrinos decouple 
and the soft parameters run as in the MSSM. The neutrino Yukawa coupling remains constant in this region. 

We find an upper limit on the magnitude of the neutrino Yukawa coupling at weak scale of 

Y N < 0.05Y t . (4.6) 

For larger values of Yjv, the model is ruled out for all values of the the cut-off scale, A/gut < A < Mpi anc k, due to a 
stau LSP over this parameter range. The large neutrino Yukawa couplings drive the left-hand lepton masses down, 
which in turn lowers the mass of the lightest stau through mixing. If the neutrino Yukawa coupling is smaller than the 
constraint in Eq. (|4.6[) . then for a region of (£, A) parameter space the results of the minimal stochastic superspace 
model can be recovered, of course now with non-zero neutrino masses due to the seesaw mechanism. Figure [3] depicts 
the mass dependence of the two LSP candidates - the lightest stau and neutralino - on the stochasticity parameter 
£ at the extremes of the cutoff parameter range. For low cutoffs, the stau is the LSP by a wide margin, however as 
the cutoff scale is increased, some phenomenologically viable regions exist for large negative values of £. Since the 
slepton mass is typically lower in the seesaw model, the acceptable region of parameter space is significantly reduced 
compared to the minimal stochastic superspace model. 



V. CONCLUSION 



The extension of minimal model with stochastic supersymmetry to incorporate neutrino masses consistent with 
experimental observations can successfully be achieved through R-parity violation. We found a specific set of trilinear 
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FIG. 3: LSP and NLSP masses as a function of £ for type-1 seesaw mechanism in stochastic superspace. These figures 
demonstrate the increased mass gap between the lightest stau and neutralino compared to the minimal stochastic superspace 
model. 



R-parity violating couplings that lead to neutrino masses which satisfy the tribimaximal mixing scheme. The char- 
acteristic feature of the resulting sparticle spectra is lighter than usual stop masses, which could lead to interesting 
decay signatures at the Large Hadron Collider. However, by choosing this method to radiatively give mass to the 
neutrinos, a natural dark matter candidate is sacrificed. This is not necessarily problematic, as the introduction of 
axions and axinos can solve this issue along with the strong CP problem, which was not addressed here. The type-I 
seesaw mechanism is an alternative approach for generating neutrino masses. We have found that it is compatible 
with stochastic supersymmetry for scenarios which have Yn < 0.051t and a generally large cut-off scale. We hope 
to investigate the viability of type-II and type-Ill seesaw mechanisms for generating neutrino masses with stochastic 
superspace in future work. 
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Appendix A: Type-I Seesaw Renormalization Group Equations 



We adopt the notational conventions of [4{ . The RGEs are the same as for the MSSM with the following exceptions 
above the seesaw scale: 



11 



dY N 

dt 
dA N 

dt 

dY v 

dt 

dA v 

dt 
dB 

~dt 

drh 2 L 
dt 



dfh 2 N 

dt 

d/j^ 
dt 
dm 2 ^ 2 

dt 



-Y, 



N 



3a 2 + -01 - 3Y/y - 4Y N 
5 



3a 2 M 2 + -O1M1 + 31^/^,7 + 4Y N A N , 
5 



^a 3 + 3a 2 + ^|ai - 6Y[/ - Yb - Y)v 
3 15 



^a 3 M 3 + 3a 2 M 2 + ^a x M x + §Y V A V + Y D A D + Y N A Nl 
3 15 
3 

3a 2 M 2 + -aiM x + 3YuAu + 3Y D A D + Y L A L + Y N A N , 
5 



3a 2 M| + -a x M( - Y L (m z L + m% + rnj li + A A L ) 



Y N (ml + m 2 N + m 2 H2 + A 2 N ) 



2Y N (ml + m 2 N + m 2 H2 + A%) 



= -/r 3a 2 + -a 1 -(3Y u + 3Y D + Y L +Y N ) 
5 



3a 2 M| + -a x Ml - 3Y V (m 2 Q + m 2 , + m 2 H2 + A 2 ,) 





Y N (m L + rh N + m H2 + A N ) 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 

(A6) 

(A7) 
(A8) 

(A9) 



where t = In (Q 2 /A 2 ). 
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